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Abstract 

In this communication we propose a most general equation to study pattern formation for one- 
species population and their limit domains in systems of length L. To accomplish this we include 
non-locality in the growth and competition terms where the integral kernels are now depend on 
characteristic length parameters a and /3. Therefore, we derived a parameter space (a, /3) where 
it is possible to analyze a coexistence curve a* = a*(/3) which delimits domains for the existence 
(or not) of pattern formation in population dynamics systems. We show that this curve has an 
analogy with coexistence curve in classical thermodynamics and critical phenomena physics. We 
have successfully compared this model with experimental data for diffusion of Escherichia coli 
populations. 
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Introduction - In recent years the phenomenon of pattern formation has been intensively 
studied to describe the spatial distribution of species in population dynamics. This amazing 
behavior of populations, observed under certain conditions, can be modeled by nonlinear 
equations of reaction and diffusion- type [l-7|. Such mathematical models provide a rich 
structure to include a variety of intra- or inter-specific interactions among species jsl, ol, as 
well as permitting to describe many forms of effects of dispersal with or without memory 



effects 
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-ll2|. 



On the other hand, the overwhelming majority of the studies have shown that Fickian- 



type diffusion 



13Hl5| is unable to describe spreading of species in population equations 



formulated through reaction-diffusion models. Moreover, this is not the only possible criti- 
cism that one can find at the ordinary nonlinear reaction-diffusion approach. In population 
dynamics context, one believes that there is no real justification for assuming that interac- 
tion among species are, in fact, local. There are many models in which such an assumptions 
become clearly unwarrantable, as for example, in competition of one-species in a habitat 
where the system is rapidly equilibrated or in typical biological interactions where the indi- 
viduals intercommunicate via chemical means. Of course, these are typical nonlocal effects 
which should not be overlooked. Other forms of nonlocal growth and interaction effects 
have been also observed when we deal with wide variety of biological fields, such as epidemic 



spread in network 



16Nl8|. embryological development, and bacterial growth (l9|, where the 



density of individuals involved are not small and the analysis of local or short-range diffusive 
flux is not sufficient accurate to understand the dynamical aspects of these phenomena. In 
these cases we need to include the contribution of long-range effects, and then to analyze 
the domains limits for the existence of patterns in these systems. 

The present communication is an attempt to build a most general equation to study 
effects of pattern formation in one-species population dynamics. Our starting point is to 
write an equation which includes nonlocal growth and interaction terms involving long-range 
effects in the system. This equation can be written as 

duix^ tl/" III f III 

= a ga{x — x')u{x' ,t)dx' — bu{x,t) / fjj^x — x')u{x' ,t)dx' , (1) 



where u{x, t) describes the population density with growth a and competition b terms. Then 
we denote by ga{x — x') the correlation growth function, which weights the growth of a 
population in the domain Q for a specific growth length parameter a. We call /^(x — x') the 
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correlation competition function which weights the interaction among the constituents of the 
population for a competition length parameter /3 in the domain Q. In the intuit to modeling 
population dynamics, we have assumed that the kernels are symmetric functions, such that 
//3 — !■ and — )■ as |x — x'| — > oo. An of advantages of Eq. ([T]) is that it provides an 
useful concept to describe a great variety of long-range diffusive effects in physical systems, 
by permitting we can absorb all higher derivatives (diffusion and dissipation terms) in only 
integral term. Moreover this model is then parametrized mainly by a and /3 quantities, 
and this allows us to suppose on the existence of values {a, /3) for which there are pattern 
formation. Indeed the relation between these length domains offers a simplest model to 
deal with quantitative estimates of experimental data related to the growth dynamics of 
bacteria, for example. In this case, specific domains which show the existence (or not) of 
patter formation, which incorporates long-range effects as well, is important for the physical 
description of spreading of individuals with nonlocal dynamics. 

Starting from Eq. ([T]) we can derive important connections with classical population 
models by carrying out appropriated limits. For example, if fp^x—x') = ga{x-~x') = 6{x—x') 
we get the logistic equation |15| 

^^^^^^^ = au{x,t) -bu{x,t)^. (2) 

One can also consider ga{x — x') with a finite range, such that we can expand the growth 

term as 

a / ga{x - x')u{x', t)dx' = 7^^r^^(^' ^) ' (3) 

where y = x — x' and the k-moments are 

¥= j y''ga{y)dy. (4) 

Using the above procedure with /^(y) = 6{y) and retaining the first two terms in Eq. ([3]) 
we get the ordinary Fisher equation 

du(x,t) ^d'^u(xA) / X , 9/ X /^s 
dt ^ dx^ ^ ''''^^'^^ " bu\x,t). (5) 

Here we show a very important point: the first gain with the nonlocal growth term is the 
possibility to connect the growth rate a with the diffusion constant 

fl = f. (6) 
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Note that this equation shows that a species with a large growth rate has "more need" 
for diffuse behavior, i.e. a large growth rate creates a large pressure proportional to the 
concentration gradient D-^u{x,t) which increases the diffusion. It shows that the diffusion 
is intrinsically related to the existence of the nonlocal growth term Qc, i.e. it is necessary to 
exist a second moment 7^ in according to Eq. (jl]). 

The higher order terms (m > 1) in expansion of Eq. Q yield the dispersive terms. It 
is interesting to note that if is not even the first derivative yields the convective term 
v-^u{x,t), where we can obtain the convective velocity as v = ya. On the other hand, if we 
have an asymmetric ga it corresponds to a convective drift |21|. If we keep the expansion 
up to second order and fpiy) in the second integral, we get nonlocal Fisher equation 

du{x, t) _ ^^^^^ t) + D ^ - bu{x, t) [ f{x - x')u{x,t)dx , (7) 



dt dx'^ jQ 

which has been widely used by many authors to discuss pattern Note that is Eq. ([1]) is 
most general, incorporating all the previous equations. 

Perturbative analysis - In the study of pattern formation through a model of popu- 



ation dynamics it is usual to calculate a quantity known as the growth rate of pattern 7 
9], [lo| which leads to the pattern formation in the system. Therefore, we first shall start 
with the perturbative analysis through the function 

u{x, t) = - + e exp (^ikx + (p{k)t^ , (8) 

where a/b is the homogeneous steady state solution, constant in space and time. The term 
eexplikx) exp(0t) is a perturbation to the steady state that will grow or die out, depending 
on the values of the wavenumbers k. Substituting Eq. ([8]) into Eq. ([T]) and retaining only 
first order perturbative terms, we find a dispersion relation between the complex pattern 
growth rate and the wavenumber k, given by 

m = a(^J^,{g^iy)}-J^,{f^iy)}^+ta(^J's{ga{y)}-J's{My)}) , (9) 

where J-^d'} and J^s{'} are, respectively, the Fourier cosine and sine transform of the influence 
function ga{y) and f/3{y) (assumed to be even). Therefore, only the real part of the complex 
growth rate = -f{k)+i^{k), where -f{k) = a{^j^go,{y)cos{ky)dy- j^f^{y)cos{ky)dy-l^, 
i. e. the growth rate of pattern 7(/i;), will be important to determine whether the perturbation 
with wavenumber k will die out or will generate a pattern, for negative or positive values, 
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respectively. Now let us consider the simple case of the square interaction influence function 
given by 

My) = ^[Q{P - y)Qif3 + y)] , (10) 

where G refers to the Heaviside function and /3 is the cut-off range (0 < /3 < L, where L is 
the size of the system). If we consider similar relation for with cut-off < a < L, from 
the condition (fTOl) . 'j{k) is given by 

7(A:) = af^^-^^-l). (11) 

\ ka kp J 

Therefore we can study self-organization of the equation ([1]) considering that the system 
depends physically on the domain of the functions and In this case, pattern formation 
appears when wave numbers A;, in growth rate of pattern, obey the condition ^(k) > 0. Note 
that for a = 0, 7(/c) is larger. This will happen for lower diffusive systems D ^ 0. 

The Eq. (fTTj) is plotted in Fig. [1] for different values of the growth length a = 
(0.01, 0.03, 0.09, 0.40) with competition length parameter /3 = 0.45 flxed. In this flgure, 
we verified that when a < /3 we may have 7(/c) > 0. If a — /3, the function jlk) becomes 
negative. This behavior of 'y{k) is very important to determine if we have a large or a 
negligible amplitude of pattern. We show in Fig. [1] that pattern formation appears for 
values a < f3 with 'y{k) > 0. This behavior is also verified later with numerical results, see 
Fig. 4, and discussed through experimental values. 
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FIG. 1: The real part of growth exponent ^{k) as a function of k plotted for different values of the 
correlation length of growth a with a length interaction of individuals f3 fixed. The 
pattern formation appear for those values of k for which 7 is positive. 
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Numerical method - To solve Eq. ([T]) numerically, we applied the Operator Splitting 
Method (OSM) (20j. By this method, the operator of the differential equation is split into 
several parts, which act additively on u{x,t). If we write Eq. ([T]) as 

du{x, t) 



dt 



Tu{x, t), 



where T is the total operator, then 



with 



Tgrowu{x,t) = a ga{x — x')u{x' ,t)dx' 
Jn 



(12) 



(13) 



(14) 



Tintu{x,t) = —bu{x,t) / fp{x — x')u{x\t)dx'. 



(15) 

In the latter equations Tg^ow and Tint are nonlocal growth and nonlocal interaction operators, 
respectively. In our numerical calculations, we have used periodic boundary conditions 
u{x = 0,t) = u{x = L,t) with spatial period L. For each part of the operator, we apply a 
known difference scheme for updating the function u{x,t) from step j to step j + 1. 
In Fig. |2]we show the evolution of u{x,t). We start with a distribution of individuals 



u{x, 0) = — exp 



[X - Xo 



\2-\ 



2a^ 



(16) 



where F 



2^ 



erf(^)+erf(^ 



and we see the evolution to a state which exhibits 



pattern. We use a = 0.3, Xq = 0.5 and L = 1.0. The spacial and time increments are 
6x = 1 X 10^^ and 6t = 1 x 10^^. In all simulations we use a = b = 1.0. Several numerical 
experiments show the final state independent of the initial conditions j2l] . These simulations 
are fundamental to show the pattern formation that appears after a long time for bacterial 
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22|. Similar simulations are used to compose Fig. [3] and Fig. |H 



growth 

In Fig. |3l we show the evolution after 60000 time steps where the density has reached its 
final form. Each curve is similar to the simulations described in Fig. [2J The first column is 
for /3 = 0.07 and for up to down a = (0.009, 0.012, 0.019). The second column has /3 = 0.11 
and for up to down a = (0.012,0.020,0.030). For each /3 the lower curve represents the 
hight value of a, for which there is no more pattern formation, i.e. for a > a*, where we 
get u{x, t) = a/b. 
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FIG. 2: The typical pattern formation on density u{x,t) as a function of x and t in arbitrary 
units. The growth rate and interaction rate are a = b = 1.0. The competition length 
parameter /3 = 0.15 and the growth length parameter a = 0.009. 

In Fig. m we show the region in the space {a, f3) where pattern can exist. For each /3 the 
points represent the a* above which there is not more pattern, such as described in Fig. |3i 
The shadow area is limited by the coexistence curve 

a*(/3) = P(/3)(/3,-/3)^ (17) 

which is the best fit of the points. is a polynomial with no roots in the region < 

/3 < /3c- For f3 ^ f3c get from Eq. (fTU]) f3c = 1/2, i.e the function fpi^y) weights equally 
the all the space < x < L, and consequently we have no pattern formation for /3 > Pc- 
Finally, from the data we have estimated the value /j, = 0.53 ± 0.06 at the vicinity of Pc for 
the exponent in Eq. (fT7|) . 

Experimental data - Starting from Eq. flTUl) we can compute = a/3. By inserting 
this result into Eq. we get 

-y?^ (IS, 

Now using the experimental values for a = (2.23 ± 0.2) x lO^^'s"^ and D = (2.2 ± 0.2) x 
10~^cm^s~"'^ obtained by N. Perry [22], for systems with Escherichia coli populations, we 
can estimate the value of a, given by a = (7.70 ±0.09)mm. Then to form pattern in a finite 
system of length L, /3 must be inside the shadow area of Fig. 4. In fact, the value of a as 
obtained in Eq. (fT8|) is determined by the coefficient of diffusion D of system and it establishes 
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FIG. 3: Snapshots of the stationary state u{x) for some values of competition length parameter f3 
and growth length parameter a with correlation competition function and correlation 
growth function Eq. (jlOp . In this snapshots we consider the growth rate a and 
competition rate b as 1.0. For a fixed /3 as q increases the pattern disappears. 




FIG. 4: The phase diagram of critical correlation growth length Oc as a function of critical 
competition length parameter f3c- The Pattern and Non Pattern region indicate the 
separation of large-amplitude patterns and negligible-amplitude patterns for Eq. (I16p . 



a lower referential limit for the presence (or not) of pattern formation. For /3 <^ L, we are in 
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the linear part of Eq. (fT71) and we get a = 0.3/3, therefore we have pattern for /3 > 25.4mm. 
In this case, only experimental values of width of influence function /3 > 25.4mm permit 
patterns, which are in concordance with our theoretical and numerical results. Moreover 
our formulation allow us to analyze pattern formation as an interplay between two length 
parameters a and /3. It is important to note that the fact that < a < /3 is not just a 
curiosity of the theory, it is one of its major result. Without a finite value of a, there will 
be no diffusion, which is fundamental for any species, and so reproduction and propagation 
are associated. Consequently one should expect a non null a. On the other hand if a is too 
large, bonds are tight, and they may face difficult to meet and consequently to reproduce. 
This phenomena described here for pattern formation in bacterial colony can be observed in 
large animals with migrate habits, such as deer and woolf, they travel with the only family, 
we call this phenomena the faithful sailor travel. 

Conclusion- The presence of memory, non-locality in time, have been used to explain 



ergodicity violation in particle diffusion 



23 



24| . Since pattern formation implies in ergodicity 



breaking, one could expect that a nonlocal space kernel would yield that. Consequently, we 
proposed here a new formulation for population dynamics, which includes a growth and 
a competitive nonlocal terms. The presence of two kernels ga{x) and /^(x) demand the 
existence of a growth length parameter a and of a competition parameter /3. Particular 
values for the kernels yield most of the previews formulations of population dynamics. We 
obtain a domain region < a < f3 where patterns may arise, a coexistence curve similar 
to those in phase transition, and a direct connection between the diffusion constant D the 
growth rate a and the mean square deviation = J fa{y)y^dy which is a function of a. 
More results can be obtained from this formulation, however there are some restrictions, we 
need more detailed dynamical experiments in growth, in such way that we can propose more 
elaborated kernels. Those present here, Eq (4), gives us a rough idea of the dynamics. More 
accurate ga{x) and //^(x) will permit us to get a better description and a generalization for 
two dimensions. 
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